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Abstract 

Single-, two-, and three-photon transitions were driven amongst five quantum states of a nio- 
bium persistent-current qubit. A multi-level energy-band diagram was extracted using microwave 
spectroscopy, and avoided crossings were directly measured between the third and fourth excited 
states. The energy relaxation times between states connected by single-photon and multi-photon 
transitions were approximately 30 — 100 /is. Three-photon coherent oscillations were observed 
between the qubit ground and fourth excited states with a decoherence time of approximately 50 
ns. 

PACS numbers: 03.67.Lx,03.65.Yz,85.25.Cp,85.25.Dq 
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Superconducting Josephson junctions are devices that exhibit quantum phenomena 
amongst their macroscopic degrees of freedom . Demonstrations of macroscopic quan- 
tum tunneling and energy level quantizationprovide a basis for utilizing superconducting 
devices in quantum computing applications |3j. Quantum-state superposition , multi- 
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have been demonstrated with Josephson-based quantum 



Qubit spectroscopy and coherent oscillations can be externally driven in several regimes: 
by single-photon or multi-photon transitions, with weak or strong driving amplitude, and 



in two-level or multi-level systems 
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single- and multi-photon quantum coherence in a strongly-driven, two-level, charge-qubit 
system [19]. Likewise, Claudon et al. observed single-photon coherent oscillations in a 
multi-level DC SQUID potential well jsjj. Recently, Saito et al. reported single- and multi- 
photon spectroscopy in a weakly-driven two- level persistent-current (PC) qubit [2]]]. Most 
demonstrations have utilized aluminum-based devices. 

In this paper, we investigate multi-photon, multi-level spectroscopy and dynamics in a 
weakly-driven niobium persistent-current qubit. We utilize static and time-dependent spec- 
troscopic information from single-, two-, and three-photon transitions to plot a multi- level 
energy band diagram of the qubit, to identify avoided crossings between its macroscopic 
states, to measure the energy relaxation time Ti from selected excited states, and to demon- 
strate three-photon coherent oscillations between the qubit's ground state and fourth excited 
state. 

A PC qubit is a superconducting loop interrupted by three under-damped Josephson 
tunnel junctions (JTJs) [see Fig. HJ^a)]. Two JTJs are designed to have the same critical 
current J c ; the third is al c . For 0.5 < a < 1 and with an externally applied magnetic flux 
fq ~ $0/2 ($0 = h/2e is the flux quantum), the system is analogous to a particle in a two- 
dimensional double-well potential with a multi-level energy band diagram as simulated (solid 
lines) in Fig. E)[b). Throughout the paper, the flux will be parameterized by its detuned 
valued 5f q = / q — $o/2- A one- dimensional slice of the double-well potential biased at a flux 
detuning 5f q ~ 3 m$ is shown in Fig.^b). Energy levels with positive and negative slopes 
correspond to macroscopic persistent currents i q of opposing sign, each associated with one 
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FIG. 1: (a) Schematic of the PC qubit surrounded by a readout dc SQUID. "X" denotes a JTJ. 
The qubit is biased by a magnetic flux f q , and the circulating current of the qubit is I q . (b) 
Simulated double-well potential of the qubit with eigenenergies Eq — E4 shown. The energy level 
separation is hv. (c) Time sequence of the microwave pulse, SQUID bias current pulse Is, and 
SQUID voltage response U3 for one measurement trial. The microwave pulse duration is r p and 
the delay time between microwave pulse off and measurement pulse on is r^. 



of the potential wells |22j. The single- well states are coupled via the potential barrier; the 
aggregate system has eigenstates with eigenenergies shown in Figs, ^b) and|2^b). Varying 
the flux Sf q tilts the double-well potential and, thereby, adjusts its eigenstates and energy- 
levels. Microwaves with frequency matching the energy level spacing can generate transitions 
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The Nb PC qubits were fabricated with a Nb trilayer process at MIT Lincoln Labora- 
tory [23J. A schematic of the qubit and readout dc SQUID circuit is shown in Fig. ^a). 
The inner loop with three JTJs is the PC qubit with a circulating current I p ~ 0.45 /^A and 
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a « 0.84. The critical current density is J c « 160 A/cm 2 , ps 300 GHz, and E c ~ 0.65 
GHz. The qubit's loop area is 16 x 16 /im 2 , with a self- inductance of L q « 30 pH. The 
readout SQUID surrounds the qubit and consists of two JTJs with equal critical current 
I c o = 2 fiA. Both JTJs are shunted with a 1-pF capacitor to lower the SQUID resonant 
frequency. The SQUID's loop area is 20 x 20 /im 2 , with a self-inductance of L$q ~ 60 pH. 
The mutual inductance between the qubit and dc SQUID is M ~ 25 pH. 

Measurements were performed in a dilution refrigerator. The device was magnetically 
shielded by four cryoperm-10 cylinders and a superconducting can. All electrical leads 
were carefully attenuated and filtered to minimize noise. For each measurement cycle, we 
initialized the system by waiting a sufficiently long time, typically 5 ms. Then, as illustrated 
in Fig. [1(c), a microwave pulse with duration r p is applied to drive quantum-state transitions. 
After delaying a time interval ra, a readout current pulse Is is sent to the SQUID, and the 
voltage Vs across the SQUID is monitored with a universal counter. The Is pulse consists 
of a 20 ns short pulse, with 5 ns rising and falling edges, and an 18 /xs trailing plateau. The 
SQUID switches to the finite- voltage state depending on the qubit state and the current Is- 
This procedure is repeated approximately 1000 times, and the SQUID-switching counts are 
used to estimate the switching probability P sw . In the following experiments, P sw = (1) 
corresponds to states with positively (negatively) sloped energy bands. 

Fig. [2(a) shows an example of how P sw varied as a function of 8f q with microwave ir- 
radiation at v — 9.8 GHz and of duration r p = 1 {is (solid black line). For reference, the 
expectation value of the qubit initial state in the absence of microwaves is fitted with the 
thermal distribution function (red dashed line) (2^. In the presence of microwaves, numer- 
ous resonances (both peaks and dips) are observed. The resonances are labelled I, Ha, lib, 
and III at the top of Fig. [2(a), corresponding to the number of photons involved in each 
transition. The flux axes in Figs. Ufa) and 12(b) are aligned to coordinate the resonances 
in the switching probability trace with their corresponding position in the energy band di- 
agram. The experimental data points in Fig. E(b) are indicated by individual markers, the 
solid lines are the energy bands obtained from simulating the qubit using the fabrication 
parameters listed above [22j, and color identifies the first six eigenenergies E . . . E 5 . The ar- 
rows in Fig. |2(b) are each of height v = 9.9 GHz; the number of stacked arrows corresponds 
to the number of photons involved in a given transition. 

The single-photon resonances, label I in Fig. 12(a), are transitions between the two lowest 
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FIG. 2: (Color) (a) The switching probability P sw vs. <5/ q for a microwave pulse at frequency 
v = 9.9 GHz with r p = 1 ^s (solid black line), and thermal distribution in the absence of microwaves 
(red dashed line). Labels I, Ha, lib, and III indicate single-, two-, and three-photon resonances, (b) 
Measured (markers) and fitted (lines) multi-level energy structure of the PC qubit. Eigenenergies 
Eq . . . E5 are shown at right. A, B, and C indicate levels represented by the three-level model (see 
text). Inset: zoom-in of the avoided crossing between levels £3 and E±. 

energy states due to the applied microwave magnetic field at frequency (Ei — E )/h = 
v. For the same microwave frequency and amplitude, two-photon resonances, label Ha in 
Fig. Efa), were also identified between the two lowest energy states at a flux bias yielding 
(Ei — Eo)/h = 2z/. In addition, a 3-photon resonance, label III in Fig. Efa), was identified 
between the ground and fourth excited state at a flux bias yielding (E^ — E )/h = 3u. 
Sweeping the microwave frequency varied the flux-distance between resonance pairs with a 
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slope |Ai//A/ q | ~ 1/n GHz/m$ , allowing us to distinguish n = 1 ... 3 photon transitions 
and, thereby, reconstruct the energy band diagram. The reconstructed energy band diagram 
is consistent with that simulated using device fabrication parameters [solid lines, Fig. |2£b)], 
and it provides clear picture of the multi- level energy-band structure of the qubit. 

The additional two-photon resonances, label lib in Fig. 12(a), are transitions between 
bands E± and Ei- A non-zero residual population may exist in E\ due to finite temperature, 
and so the lib-type transitions were generally observed for flux biases at which E\ had non- 
zero equilibrium population [dashed red line in Fig. |2t a )]- Since E 2 has the same slope and, 
thus, circulating current as the ground state energy Eq at these flux biases, P sw decreased 
(increased) for Sf q < (8f q > 0). This is in contrast to transitions I, Ha, and III, for which 
P sw increased (decreased), because the accessed bands E 1 and E± have opposing slope and 
circulating current to Eq. 

At avoided crossings, the energy bands have zero slope, and the persistent current van- 
ishes. Thus, at 5f q = 0, the SQUID switching- current readout cannot distinguish the Eq 
and Ei states. However, we estimated the Eq and Ei tunnel splitting Aj.2 < 0.1 GHz by 
fitting to the eigenenergy difference E\ — Eq = \/Af 2 + e 2 , where e is the flux-dependent 
energy difference between the uncoupled states [5J]. The simulated energy band diagram 
yielded A 12 ~ 10 MHz. The MHz tunneling splitting indicates that the coupling between 
the £0 and Ei states is very weak. 

We directly observed an avoided crossing with a tunnel splitting A34 « 500 MHz be- 
tween energy levels E3 and E± using multi-photon transitions [inset Fig. Efb)], indicating a 
superposition between two macroscopic quantum states. Although the slopes of £3 and E± 
are zero at the avoided crossing, their SQUID readout signals are distinct from that of the 
ground state, which has finite slope at this flux bias. By driving transitions from the ground 
state, we could directly measure the E 3 — E4 avoided crossing and its tunneling splitting 
A 34 . 

The energy relaxation times T\ between macroscopic quantum states were measured by 
driving population to the excited states (E\ and £4 in Fig. EJ) and monitoring it as a function 
of Td- Fig.|3fa) shows examples of the population decay; an exponential function fits the data 
remarkably well. The obtained relaxation time from Ei is Ti ~ 30 /is, and the relaxation time 
from £4 is even longer, T\ ~ 80 fis. Both are consistent with the independently measured 
intra- well relaxation time for a similar qubit design j^J, 13, |^. We also investigated 7\ 
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FIG. 3: P sw at v = 9.9 GHz as a function of the readout delay time t^. Using a fit exponential 
decay, we obtained T\ ~ 30/is for the I transition, and T\ ~ 80 /is for the III transition in Fig. [21 
Inset: T\ as a function of v for the I transition. 

as a function of the energy difference hv (inset Fig. |3J). The experimental 7\ tends to 
increase with increasing z/, qualitatively agreeing with the spin-boson bath model Q, ll^ . 
However, the deviation between the data and theory suggests that we may have other 
sources and structure to the dissipation. Nevertheless, we can potentially benefit from this 
long T\ in future e„ ts , such a, the three-leve! te qublt and superconductive 
electromagnetic induced transparency 

Although the coupling between the ground and first excited states proved too weak to 
see definitive Rabi oscillations, we were able to observe oscillations in P sw for transitions 
between the ground and fourth excited states, since the presence of the avoided crossing 
(Fig. El inset) gives a larger overlap between these states. We biased 5f q at the 3-photon 
resonant peak (III in Fig. |2J) and measured the population in the excited state P sw as a 
function of microwave pulse duration time r p , shown in Fig. Ilfa)-(c). The population P sw 
oscillates as a function of r p , and the oscillation frequency changes with the microwave 
amplitude. The Rabi decay time is approximately 50 ns, although the oscillation visibilities 
are only approximately 20%. In Fig. 4(f), the oscillation frequencies were plotted as a 
function of the microwave amplitude. Coherent oscillations were also observed at several 



7 




Pulse Duration t p (ns) Amplitude (a.u.) 

FIG. 4: (a)-(c) The amplitude of the 3-photon resonant peaks (III in Fig. EJ) vs. the microwave pulse 
duration r p for v = 9.90779 GHz at different powers. The population oscillations with different 
frequencies are clearly observed, (d) Coherent oscillations simulated using a three-level model, 
(e) Lorentzian-fitted (red line) resonance peak data (blue markers) with FWHM ~ 5 MHz. The 
corresponding dephasing time is ~ 100 ns. (f) The Rabi frequency as a function of the microwave 
amplitude. The solid (dashed) line is a best fit to a third-order (first-order) Bessel function of the 
first kind. 

frequencies in the range 9 < v < 10 GHz near the avoided crossing (not shown). 

The oscillations in Fig. 0] arise from multi-photon transitions in a multi-level system. 
We modelled the multi-level, multi-photon transition using a three-level system in which 
coherent oscillations are driven between levels A and B, and level B is taken to be strongly 
coupled to a third level C [see Fig. El^b)]. Fig. IHd) shows coherent oscillations simulated 
using this model, assuming coupling strengths and decoherence parameters consistent with 
our qubit. Details of this model will be presented elsewhere, and we only state the qualitative 
results here: driving levels A and B modifies the coupling between B and C, and this 
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modification is dependent upon the driving amplitude. The presence of the third level 
C results in an amplitude-dependent visibility that vanishes for small driving amplitudes, 
while, for larger amplitudes, the power-dependent coupling to the third level C effectively 
detunes the A - B transition. The net result is a limited window of driving amplitudes over 
which multi-photon transitions can be clearly observed in a multi-level system. 

In the limit that the coupling to C is eliminated, the above model reduces to the n- 
photon, two-level case, in which the Rabi frequency will vary as J n {x), where J n is the n-th 
order Bessel function of the first kind and x = eV^jhv is proportional to the microwave 
driving amplitude V ac [19] . The experimental amplitude dependence of the Rabi frequency 
in Fig. BJf) is reasonably fit by a third-order Bessel function over the limited range of ampli- 
tudes, with the reduced visibilities observed experimentally in Fig. Efa)-(c) and simulated 
in Fig. Efd). Despite the fringe-contrast reduction due to the multiple levels, the ~ 50 ns 
decoherence time of our Nb PC qubit is comparable to those of the qubits fabricated with 



In summary, we demonstrated multi-level, multi-photon coherent dynamics in a niobium 
PC qubit. From the single- and multi-photon resonances, we extracted a multi-level energy 
band diagram. We directly measured an avoided crossing, and we demonstrated three- 
photon coherent oscillations between the ground and fourth excited state in a multi-level 
environment. 
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